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Abstract 

A method is given to obtain closed form formulas for the energy and forces for an aggregate of 
charges interacting via a logarithmic interaction under periodic boundary conditions. The work 
done here is a generalization of Glasser's results [M. L. Glasser, J. Math. Phys. 15, 188 (1974)] 
and is obtained with a different and simpler method than that by Stremler [M. A. Stremler, J. 
Math. Phys. 45, 3584 (2004)]. The simplicity of the formulas derived here makes them extremely 
convenient in a computer simulation. 
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I. INTRODUCTION 



Numerical simulations are routinely employed in the study of physical problems that are 
difficult to solve analytically. Since it is not possible to simulate realistic physical systems, 
containing ions of the order of Avogadro number, one usually works with a very small system. 
For small systems, containing a few hundred to a few thousand charges, boundary effects 
become relatively pronounced, especially if the nature of interaction is long range. To avoid 
this problem, periodic boundary conditions (PBC) are usually employed. 

In many simulations, the nature of interaction is such that the potential satisfies the 
Poisson equation. For example, a logarithmic interaction in two dimensions (2D) satisfies 
the Poisson equation in 2D. The fact that these potentials satisfy the Poisson equations 
makes it easier to treat them in the Fourier space. Our aim in this paper is to consider one 
such especial case. We obtain formulas for the energy and forces for a system of particles 
interacting via the long range logarithmic potential under PBC. Related work on the log- 
arithmic potential could be found in the papers by Gr0nbeck- Jensen*, Liem-. and Tyagi^. 
However, all these works do not give formulas in a closed form. Recently, Stremler— general- 
ized Glasser's^ work to obtain formula for a kind of lattice sums that have a direct bearing 
on the logarithmic interaction under PBC. In this paper, we take a different approach to ob- 
tain Glasser's and Stremler's work and show the connection of their work to the logarithmic 
interaction under PBC. The derivation given here generalizes Glasser's work for the case of 
a rhombic cell. The method presented is much simpler and straightforward than that given 
by Stremler— . The outline for the rest of the paper is as follows. In Sec. II we present the 
method to obtain the closed form formulas, and in Sec. Ill we discuss the results. 

II. CLOSED FORM FORMULAS 

For two charges interacting with a pair-wise potential under PBC, the interaction of a 
charge with the other one consists of three parts. The first is the direct interaction between 
the charges. The second part consists of the interaction of the first charge with all of the 
periodic images of the second particle, and the last part consists of the interaction of the 
first charge with all of its own periodic images. It does not matter which of the two charges 
we call the first or the second charge. These three interactions when summed over all pair 
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of charges contained in the simulation cell form a series. There are two important points 
about this series. The first is that the series will diverge if the system is not overall charge 
neutral and the second that even if the system is charge neutral the series would still be only 
conditionally convergent. In the simulation of a physical system, one is interested in the 
interaction energy of the charges in the basic simulation cell with all the charges in a big box 
of size R x R. The box is chosen such that the basic simulation cell is located at the center 
of the box. It is clear that for a charge neutral system the result will converge to a well 
defined number as R tends to infinity. One way of obtaining this desired limit is through 
the use of artificial background charges. The idea is to consider that every single charge 
q, and its images come in conjunction with a uniform distribution of charge such that the 
total charge, in any layer, due to this uniformly charged layer adds up to —q. For a charge 
neutral periodic system, imposing these uniform background charge distributions for every 
charge in the system does not matter since the total uniform background charge adds up to 
zero at every point of the space. However, now a charge located within the basic simulation 
cell at position r may be thought of not only interacting with a second charge located at the 
origin and its periodic images, but also interacting with the neutralizing background charge 
of the second particle and its images. Thus, if we have N charges located at positions r< 
in the simulation cell, then the interaction is given by 



E toUl = E E MiG (x) + ^ E ^if, (2.1) 

i=l j=i+l i 

where = 1% — Tj, and A has arbitrarily been chosen as the length scale and G(r/X) 
satisfies the modified 2D Poisson equation 

V> G (^-f ? *(^)+^, ,, 2 , 

where l\ and l 2 denote the lengths of the sides of the rhombic cell, and < 6 < tt is the 
angle between the two unit vectors eiand e 2 , characterizing the unit rhombic cell with 
ei.e 2 = cos 6. The vector r and I are given by 

r = r 1 e 1 + r 2 e 2 , I = n\e.\ + n 2 e 2 , (2.3) 

where n\ and n 2 are integers, each ranging over —00 to +00. The second term on the 
right hand side of Eq. (j2.2j) amounts to the presence of a neutralizing background charge 



mentioned above. The G part represents the interaction between two different particles, 
while the G se if part describes the interaction of a particle with its own images. One notes 
that Ggeif may be easily obtained from G as 



G S eif = lim 
M->o 



G, D +ln ( A 



(2.4) 



The force components may be obtained from G as well by taking the gradient with respect 
to the position of the particle. Thus, the whole problem reduces to obtaining G. The first 
step in this direction is to realize that the solution of Eq. (|2.2jl in the Fourier space is given 
by 

„(\r\\ 2tt /V^ expUQ ■ r) l\ , , 

where 



Q = n 1 b 1 +n 2 b 2 . (2.6) 

Here, because of the PBC, we may assume |ri| < h/2 and \r 2 \ < Z2/2. The sum over Q runs 
over all reciprocal lattice vectors spanned by 

2tt 



(ej — e,j cos 6), (2.7) 



* li sin 2 9 ' 

for = (1, 2), (2, 1) and n\ and n 2 are integers. From now onwards we will assume that 
a final limit £ — > is to be taken. At this point we note that the modified Bessel function 
of the second kind satisfies the relation 

(v > + n Ko (M) = (x) • < 2 -«> 

This means that the solution of 
may be simply written as 

On the other hand, we can obtain the solution of Eq. ()2.9|) directly in the Fourier space by 
substituting 

Q 
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and using the identity 



^ V A J hhsmU^ 



Q 

This results in 

t\ 2n exp(«Q.r 



Thus, using Eqs. (12. 5J) and (j2.13|) we can express G (r/X) as 

G G)^ m o( G »(f-)-y^bF)' (2M) 

which by making the variables r and £ dimensionless through scaling with length A and 
using Eq. (|2.1U|) can be written as 



Hi 



A/ [ V A/ asin^ 2 

where a is ^h/X 2 - From here onwards we will choose l\ to be our length scale thus we will 
set A = l\. Now writing 

j y j ^ | ^ ^ 1/2 

— - — — [(x + m) +(y + na) + 2 (x + m) (y + no) cos 9] (2-16) 
= [{x + m + (y + n)a cos 6} 2 + (y + n) 2 a 2 sin 2 0] 1/2 , (2.17) 

and using the identity 4 

m=-oo ^ ' m=i i/a 2 + (27rm) 

x exp i^~ r \J Oi 2 + (27rm) 2 ^ cos (27rmx) , (2.18) 

we can write 

exp (— ^crsin^ |y + n|) 



i n=— c 



+ — exp (— 2-nma sin 6 \y + n\) 

7i=— co m=l 

x cos [27rm {a; + (y + ra) a cos 0}] , (2-19) 



oo 

oo oo 
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where x = rx/h, y = r 2 /h and we have substituted £ = in the second part on the rhs of 
Eq. (|2.19p . The first sum over variable n in Eq. (j2.19|) can be carried out as it is just a 
simple geometrical sum: 

^ exp(-£q sinfl \y + n\) _ a sing cosh [a(l-2|y|)] 

^ £ ~ 2 asinhta] ' 1 j 

n=— oo L J 

where a = (cr£ sing) /2. Now we note that 

/ cosh \a (1 — 2 It/DI 1 \ 1 , . . ,9 N , 

im 1 v , . , m - — = - (1 -6b + 6k/ 2 ) . 2.21 

a->o V asinh[«] a 2 J 3 v m m ; v 7 

Thus, taking the required limit in Eq. (J2.15j) with the help of Eq. (J2.21j) we can write 

^cx siii6> (1 — 6 |y| + 6 |2/| 2 ) (2.23) 



6 

00 00 



+ — exp (— 2ixmz n ) cos i2 r nmy n ) , (2.24) 



m 

n=— 00 m=l 



where 

2 n = a sin # |y + n| , (2.25) 

y n = x + (y + n) a cos (2.26) 

The sum over m in Eq. (j2.24j) can now be carried out as it just amounts to a simple 
geometrical series. In fact, we can write 7 

y - exp (—2npz) cos (2npy) = — Re [In {1 — exp (— <r)}] 
p=i 

= -In |1 -exp (-c)| for2>0, (2.27) 

where = 2ir(z±iy) with the choice of ± sign is up to us. Now, for all n > 1, while 
summing over m in Eq. (|2.24j) . we will set q n = 2tt (z n + iy n ) where 

Vn= {x + (y+\n\)cr cos 6} (2.28) 

z n = a sin 9 \y + n\ = a sin 9 (y + \n\) , (2.29) 

and for all n < 0, we set <j n = 27r (z n — iy n ) , where 

z n = a sin 6* |y + n\ = a sin 9 (—y + |n|) , (2.30) 

y n = {x + (y - \n\)cr cos 9} . (2.31) 



Combining together the terms corresponding to n and —n, we can write the last part of Eq. 
(EHH) as 



oo oo ^ 

/ (x,y) = — exp (-2nmz n ) cos (2^my n ) 

L — ' L — ' m 

n=—oo m=l 

oo 

= — In |1 - exp (-? )| - y]ln(|l - exp (-<r +n )| |1 - exp (-<T- n )|) 

71=1 
oo 

= - In 1 1 - exp (2wiz) | - ^ In 1 1 - 2q 2n cos (2vr;z) + q 4n | , 



71=1 



where 



and 



z = x + ya exp (zt 



g = exp [nai exp (?£ 
Using the definition of the Jacobi theta function 8 : 



(2.32) 

(2.33) 
(2.34) 

(2.35) 
(2.36) 

u^O sinu u->0 cos it 
where a prime over the Jacobi theta function denotes the differentiation with respect to the 
first argument. Thus, we obtain 

<(0,g) = lin / l(M ' g) 
u^o sm u 

oo 



7?i (it, g) = 2g 1/4 sin u JJ (l - 2q 2n cos 2u + g 4n ) (l - g 2n ) 



n=l 



we note that 



lim tfl , (ti ' g) =lim^ (tt ' g) =<(0, g ) 



2q 1 ^Y[(l-2q 2n + q in ) {l - q 2n ) 

71=1 
/ OO 

2? i/4 n (i - ^ 2n ) 



^2.37) 



.71=1 



So, we can write Eq. ()2.35|) as 



I [ (l -2g 2n cos2w + g 4n ) 



71=1 



2 -2/3 g -l/6 ^ g ) 



^2.38) 



Using Eqs. ()2.32|) and (|2.38|) we can write 



/ (x, y) = h 7r \y | a sin 9 na sin 9 — In 

3 6 



■&! (tcz, q) 



ft (0, ff)] 



1/3 



^2.39) 



Thus, using Eqs. flUHD , $£32$ and (jOty we finally obtain 



G 



, ,2 1 

Trasm9\y\ In 2 — In 

\y\ 3 



>& x (ttz, q) 



1/3 



(2.40) 



The result in Eq. (j2.40j) is in agreement with that obtained by Stremler—. An especial case of 
the method gives Glasser's 6 results. As noted earlier by Stremler, there is a trivial mistake 
in the paper by Glasser. 

The self-energy of this system can now be easily obtained by employing Eq. ()2.4j) . First 
we note that as z — > 0, we have i?x (ttz, t) tf'j (0, g) ttz, where we have used Eq. f!2.37|) 
along with the fact that sin (ttz) ~ ttz. Also, it is easy to see that 

\z\ — \x + yaexp (iQ)\ 

= [x 2 + y 2 a 2 + 2xya cos P 1 



\/r\ + rf + 2rir 2 cos 6 
h ' 



(2.41) 



The self-energy is thus given by 



G 



self 



lim 

n— >0,r2— »0 



G 



In 



i/r^ + r| + 2rir 2 cos 6 



— In 2 — In tt — In 

3 



<(0,g) 



2/3 



^2.42) 



The expressions for the force may be derived by differentiating G (r/X) with respect to the 
Cartesian coordinates x'and y', which are related to the rhombic coordinates r x and r 2 by 



ri — x' — y' cot 9 
f'2 = y'/ sin 6. 



(2.43) 
(2.44) 



Thus, we obtain 



F x , (r) 



d / r \ dri d / r\ dr 2 
dri \hj dx' dr 2 \hj dx' 



TT 
h 



Re 



' ■d\ (ttz,t) ' 
(vrz,r)_ 



(2.45) 



S 



and 





d ( r \ dri d ( r\ dr 2 
dri \h J dy' dr 2 \hj dy' 
7i cot 6 i?! (nz, r) 2tx 




2n 7r (irz, t) 



(2.46) 



/i h y&x{ixz,T) 



We have thus obtained complete expressions for G(r/X), the self-energy, and the forces. 
These expression for the energy, the self-energy and the forces were numerically checked and 
they show perfect agreement with our earlier results^. 

III. DISCUSSION 

The method given here has its origin in the beautiful work of Glasser- who obtained 
the results for logarithmic interaction in a closed form but did so only for a orthorhombic 
simulation cell. However, he also reminded his readers that the results could be easily 
obtained for a rhombic cell. Although, his work was directly related to the logarithmic 
interaction, it was not noticed by many researchers working on logarithmic interaction under 
PBC in recent years. The reason behind why this beautiful work went unnoticed is that 
it was published under the category of lattice sums, with the word logarithmic not even 
mentioned in the whole paper. Stremler— has recently extended Glasser's work to include 
the case where the simulation cell may be rhombic. But the derivation given by him is 
tedious. In this paper we have derived direct formulas for the energy and forces on particles 
interacting via the long range logarithmic interaction under PBC. The method presented 
here is much simpler than that given by Stremler. 

The advantages of the formulas derived here lies in their simplicity and easy implemen- 
tation in a simulation. Most functions used in the expressions for energy and forces, such 
as the Jacobi Theta function, are already implemented in most libraries of mathematical 
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interest . 
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